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Abstract. Using multiple Bernoulli series, we give a formula in 
the spirit of Euler-MacLaurin formula. We also give a wall crossing 
formula and a decomposition formula for multiple Bernoulli series. 
The study of these series is motivated by formulae of E. Witten for 
volumes of moduli spaces of flat bundles over a surface. 
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1. Introduction 

Consider a sequence of vectors $ lying in a lattice A of a real vector 
space V. We denote the dual lattice of A by T. Let r reg ($) = {76 
r| (0,7) 7^ 0, for all G $} be the set of regular elements in T. 
Let Z be the fundamental domain in V for A and dv the Lebesgue 
measure giving measure 1 to Z. Here in introduction we freely identify 

distributions and generalized functions via this choice of dv. 
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In this paper we study the distribution £>($,A) on the torus V/A 
defined via its Fourier coefficients as: 

B^,A)(v)e-^^dv = hhJ^-r) lf 7 G 
z [0 otherwise. 

We have 

„(2i7TD,7) 

B(*,A)(«) v 



This sum, if not absolutely convergent, has meaning as a distribution. 

We call £>($, A) the multiple Bernoulli series associated to $ and 
A. They are natural generalizations of Bernoulli series: for A = 
and $fc = [to, to, . . . , to], where u is repeated k times with k > 0, the 
distribution 

p 2iirnt 

is equal to — ^B(k,t — [t]) where B(k,t) denotes the k th Bernoulli 
polynomial in variable t. 

Multiple Bernoulli series appeared in the work of E. Witten in the 
special case where the sequence $ is comprised of positive coroots of a 
compact connected Lie group G with multiplicity 2g — 1 and A is the 
coroot lattice of G. Witten shows that ([17], §3) for the above instance 
of <3? and A and for a regular element v £ Z, the value of £>($, A)(t>) 
is (upto a scalar depending on G and g) the symplectic volume of the 
moduli space A4(G, g,v) of flat G-connections on Riemann surface of 
genus g with one boundary component, around which the holonomy is 
determined by v. 

For example, consider G = SU(3), denote its simple roots by {a±, a^} 
and associated coroots by {H ai , H a2 }. Then on a Riemann surface of 
genus g = 2, the symplectic volume of the moduli space of flat SU(3)- 
connections with one boundary component marked by v — aiH ai + 
a 2 H a2 (lying in the fundamental alcove) is given by the following sum 

p 2i7r(niai+n2a 2 ) 
V 

(2i7™ 1 ) 3 (2z7m 2 ) 3 (2i7r(n 1 + n 2 )) 3 

ni^O, n27^0, rti+n 2 7^0 

up to a scalar multiple. 

These series have been extensively studied by A. Szenes ( [T2] . [T3] ) . 
who gave multidimensional residue formulae for them. 
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If / is a function on the real line, smooth and sufficiently decreasing, 
also with sufficiently decreasing derivatives, then the Euler-MacLaurin 
formula gives 



where denotes the k th derivative of /. 

We give a natural generalization of this formula involving A) in 
Theorem 16.11 The difference between the discrete sum ^] AeA /(A) and 
the integral of / over V will only involve derivatives (FLey^)/ over 
'long subsets' Y of $, that is, subsets Y such that their complement in 
$ do not span the vector space V. 

We start with giving some properties of the distribution B{<&, A) that 
are pertinent for what follows. 

The distribution $($, A) is periodic with respect to A. Moreover, it 
satisfies a certain recurrence relation which we outline next. 

We will call a set with multiplicities a list. Suppose $ contains with 
multiplicity m, then we denote the list that contains with multiplicity 
m — 1 by $ — {0}; whereas we denote the list where all copies of are 
removed by $ \ {0}. More generally, for a subset B of V, by $ \ B we 
mean the list of elements of $ not lying in B. By $ fl B we mean the 
list of elements of $ lying in B. 

For an element in $, we associate two lists of vectors as follows: 
First, we consider the list $—{0} in V, and respectively the distribution 
£($ - {0},A) on V. Let V := ^/M0 and let A denote the image of 
the lattice A in Vq. Secondly, we consider the list $o of elements of 
Vo consisting of the images of the elements in $ — {0}. Then we may 
consider B(<&o, Aq) as a distribution on V 'constant in the direction 
of 0'. Observe that if $ contains more than one copy of then $o 
contains the zero vector and consequently £>($o, Ao) is identically zero. 

It is clear that the distribution A) satisfies the following recur- 
rence relation 



Assuming that $ = [0i, 2 , . . . , 4>n] spans a pointed cone, we may 
define the tempered distribution T($) defined on test functions / by: 




(1.0.1) 



fyB($, A) = £($ - {0}, A) - 8($o, A ). 



(1.0.2) 
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Figure 1. T([e x , e 2 , e x + e 2 ], Ze x © Ze 2 ) 

In other words, T( [</>]) is the Heaviside distribution (H(<fi), f) = / °° f(t<p)dt 
and T($) is the convolution product 

H(4>i)*H((f) 2 )*---*H(4> N ) 

of the Heaviside distributions H(<pk). 

If $ generates V, T($) is a positive measure on V given by integra- 
tion against a piecewise polynomial function called a multispline. For 
any </>€$, 

(1.0.3) ^T($)=T($-{0}). 

We remark the similarity between the recurrence relations fll.O.lf) and 
(jl.0.3p . In fact we will express B(<b,A) in terms of superposition of 
multispline functions in Theorem 19.31 

If $ generates V, then the periodic function £>($, A) is piecewise 
polynomial; this we reprove in Section [51 

In the rest of the introduction, for simplicity, we assume that $ 
generates V. We call a connected component of the complement of 
affine walls (that is, hyperplanes that are generated by some elements 
of $ and their translates with respect to the lattice A) in V a tope. 
For example, Figure [1] depicts topes associated to the system $ = 
[ei, e 2 , ei + e 2 ] and A = Zei © Ze 2 . 

Given a tope r associated to the system ($, A), we denote by Ber($, A, r) 
the polynomial function on V such that the restriction of B(<&, A) to r 
coincides with the restriction of Ber($, A, r)(w) to r. 

Let W be an hyperplane in V spanned by some elements of and 
let E G T be an equation of this hyperplane. We reverse the directions 
of 'half of the 0, in $ \ W in order that they all lie in the strict half 
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space determined by E, and define 

T($\W,E):= 11 -Hi-fa)* J] Hifa). 

(<f>i,E)<0 (<t>i,E)>Q 

T($ \ W, E) is a distribution supported on E > 0. We similarly 
define T($\W,-E). 

Now we compare the polynomials Ber($, A, r) associated to two ad- 
jacent topes separated by an hyperplane W (cf. Section [7]). The jump 
can be expressed in terms of a lower dimensional multiple Bernoulli se- 
ries and a multispline function. More explicitly, we have the following 
wall crossing formula: 

Theorem 1.1. Let T\ and r 2 be two adjacent topes separated by the 
hyperplane W defined by E with (v, E) > for any v G t%. Denote by 
T\i the tope with respect to the system ($r\W,AC\W) containing tTHtJ 
in its closure. Let Ber T12 := Ber($flW / , AflW 7 , Tyi)dh be the polynomial 
density on W determined by tvi- Then, 

(Ber($, A, n)-Ber($, A, r 2 ))dv = Ber Tl2 *T($\W, £)-Ber Tl2 *T($\W, - 

The left hand side of the above equation is a polynomial density; it 
is easily proven that the right hand side is also a polynomial density. 

The wall crossing formula given in the above theorem is analogous 
to the formula in Boysal- Vergne p] . This formula is also similar to wall 
crossing formulae in Hamiltonian geometry for the push-forward of the 
Liouville measure; indeed, when crossing a wall, this piecewise poly- 
nomial measure changes according to the same scheme [H], [7J- Our 
wall crossing formula in Theorem 11.11 is thus in accordance with the 
fact that for special cases £>($, A) computes the volume of the moduli 
spaces Ai(G, g, v). These spaces can be described as symplectic reduc- 
tion at v of the Jeffrey-Kirwan extended moduli space M(G, g), so that 
their volume at v is given by the push-forward of the Liouville measure 
on M(G,g), a piecewise polynomial function periodic with respect to 
a lattice A. Recall that Jeffrey-Kirwan ([8]) proved wall crossing for- 
mulae for integrals on moduli spaces Ai(G, g,v), and used them in a 
fundamental way to compute intersection pairings on A4(G, g, v) when 
G = SU(n). However, in the general situation that we are considering 
here, we do not have such a geometric interpretation of the multiple 
Bernoulli series. 

In Section [S] we generalize the above results to the case of affine 
arrangements. 

In Section [9], we give a decomposition formula for £>($,A), describ- 
ing it as a superposition of 'basic pieces' made of convolution products 
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(a) (b) 
Figure 2. The decomposition of <B($ 2 , A)(tw). 

of lower dimensional Bernoulli series and multisplines. More precisely, 
we say that s is an admissible subspace of V if s is spanned by some 
elements of and we say that a is affine admissible, if a is a translate 
by A of an admissible subspace. Given a tope r, we express the dif- 
ference between the piecewise polynomial density £>($, A) and the 
polynomial density Ber($, A, r) as a sum of distributions A, a, /3) 
associated to proper affine admissible subspaces a and the choice of an 
element /3 G r. The supports of these distributions do not intersect 
t and are convolution products of polynomial distributions supported 
on a with multisplines distributions directed towards the exterior of 
r. Our construction is inspired by the stratification of a Hamiltonian 
manifold M using the square of the moment map as Morse function, 
and we will use a scalar product on V. Our decomposition formula is 
very similar to Paradan's decomposition of the equivariant index of a 
twisted Dirac operator on M [IT] . In [15], Paradan's decomposition 
was proved by combinatorial methods, and used to give a proof that 
quantization commutes with reduction for compact Hamiltonian mani- 
folds. We follow here very closely the line of approach of [15] . However 
our superposition is an infinite (but locally finite) superposition. This 
is in accordance with the fact that for some special cases, our distri- 
butions are related to Liouville measures of noncompact Hamiltonian 
manifolds such as M(G, g) with infinite number of critical components 
for the square of the moment map. 

For example, the periodic polynomial —B(2,t— [t])/2 in Figure Et^a) 
is decomposed in Figure El^b) as a superposition of a polynomial density 
with an infinite number of spline functions. 

Acknowledgements. We thank Michel Duflo for various suggestions. 
The first author wishes to thank Ecole Polytechnique (Palaiseau, France), 
and the second author wishes to thank Bogazigi Universitesi (Turkey) 
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TJ 


f np nnpil of TX" i* TJ 


(, > 


the pairing between [/ and V . 


r 


a lattice in U; 7 G T. 


A := r* 


dual lattice in V; (T,A) C Z, A G A. 


$ 


a sequence of vectors in Vq; G 


B(*,A) 


multiple Bernoulli series associated to $ and A (equation 12.0. 3D. 


B(k,t) 


k th Bernoulli polynomial. 


T($,A) 


the set of topes associated to the system ($, A); r a tope. 


T(X) 


multivariate spline distribution defined for a set of vectors X in V. 


n 


set of ^-admissible subspaces of V. 


w 


^-admissible hyperplane. 


H<f> 


hyperplane in U comprising of vectors u satisfying (u, (f)) = 0. 


U = U{<&) 


hyperplane arrangement associated to $. 


n H 


ring of rational functions on U with poles along %. 


Qu 


a subspace of H<n defined in |5.3[ 



2. Multiple Bernoulli series and hyperplane 
arrangements 

Let Vq be an r-dimensional vector space over Q, and let V be the 
real vector space Vq R. Let U denote the dual vector space to V. 
Let A be a lattice in V contained in Vq and r C U be the dual lattice 
to A so that (r, A) C Z. For a subset S of V, we denote by < S > the 
subspace of V generated by S. 

Let Z denote the fundamental domain in V for A. Let d^v be the 
Lebesgue measure on V giving measure 1 to Z. Our main object of 
study is certain piecewise polynomial densities on V. For our purposes 
it will be convenient to use the language of distributions. If / is a 
locally L 1 function on V, or more generally a generalized function, 
then f(v)d\v is a distribution on V. We use the notation f(v)dAV, 
although the value of / at the point v G V has usually no meaning. 

For t>o G V, the translation t(t>o) acts on distributions on V. If 
D = f(v)d&v, then t(v )D = f(v + Vojd^v. We identify a distribution 
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D on V periodic with respect to A (that is t{X)D = D for any A G A) 
to a distribution D on the torus V/A. 

We will say that a locally L 1 function / is piecewise polynomial, if 
there exists a decomposition of V in a union of polyhedral pieces 
such that the restriction of / to Cj is given by a polynomial formula. 
We then say that the distribution f{v)d^v is piecewise polynomial. 

If v G V, we denote by 5 V the 5 distribution at v: (5 V , f) = f(v). 
The Poisson formula reads as the following equality of distributions 

(2.0.1) j2 s * = J2 e2inb,v)dAV - 

We now introduce the main object of study of this article. 

Let $ be a sequence of vectors in Vq. Let 

£/re g ($) = {ueU\ (<j),u) ^ 0, for all (j) G $}. 

We will denote T fl U TCg (&) by r rcg ($). Consider the distribution 
i3($,A) on V/A defined via its Fourier coefficients: 



(2.0.2) / JB($, A){v)e-^> = { ™ ^ 7 G 

otherwise. 



We then have 



e (2i7TO,7) 



(2.0.3) £($,A)= V = ; -^d A v. 

7 er fe g(*) 11 *e* x r ' 

The above sum, if not absolutely convergent, is defined as a distri- 
bution. We call A) the multiple Bernoulli series associated to $ 
and A. Clearly, it does not depend on the order of the elements <fi in 
the sequence $ (it only depends on $ multiset). 

Remark 2.1. The formula (12.0.21) for the Fourier coefficients of B($>, A) 
is very similar to the formula for the Fourier transform of the multi- 
spline distribution T($) (defined in (jl.0.2p ) on V: if $ spans a pointed 
cone, then the Fourier transform of T($) is a generalized function on 
U satisfying 

T($)(v)e- {2i7TV > x) ' 



]l0 G $ (2ivr0,x) 



on the open set of U given by n</>G* (^j ^) 7^ 0. 
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We now list some properties of the distribution £?($, A): 

• If $ is the empty set, then 

7£r AgA 

is the ^-distribution of the lattice A. 

• If $ contains the zero vector, then A) is identically equal to 
zero. 

• A) is periodic with respect to A. 

• Let Ai C A 2 . Then i3($,A 2 ) is obtained from Ai) by averag- 
ing over A 2 /Ai. 

(2.1.1) £($,A 2 )= Yl t(A 2 )S($,A 1 ). 

A 2 eA 2 /Ai 

The above relation follows from the fact that, if 7 G Ti \ T 2 , then 

Z^AaeAa/Ai e ~ U ' 

• The distribution £>($, A) is supported on < $ > +A. 

Indeed, it is immediate to verify that (1 — e( 2i7 ™'^)i3($, A) = for 

all v g rn < $ >- L . 

• If $ generates V^, then £>($,A) is piecewise polynomial. We will 
give a proof of this property in Section 

When the lattice A is fixed, we often use the measure d^v to identify 
distributions and generalized functions. 

Example 2.2. Let A = Zw, and let = [u,u,...,u] where u is 
repeated k times. If k = 0, then #($ fe , A)(tu) = Y. n ez e2inntdt is tlie 
5-distribution of the lattice A by Poisson formula. If k > 0, then 

B{* k ,A)(tu>) = ^——dt = -jjB(k,t-[t])dt, 

where B(k,t) denotes the k th Bernoulli polynomial in variable t and [t] 
is the integer part of t. (Our normalization for the Bernoulli polynomial 
is that of Maple). In particular, for k — 1, we have B($i, A)(tu)dt = 
(\-t+[t])dt (see Figure [3]). 

If k > 1, the above series is absolutely convergent and £>($fc, A)(£a;) 
is given by integration against a continuous function on M. 

Example 2.3. Let V = Re x © Me 2 with lattice A = Ze x © Ze 2 . Let 
$ = [ei, e 2 , ei + e 2 ]. We write t> G V as t> = t^ei + f 2 e 2 . 
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\ 2 




-i/h 
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\ 2 \ 3 



Figure 3. Graph of the function B($i, A)(tu) = (§ - t + [£]) 



Figure 4. T{[e x , e 2 , e x + e 2 ], Ze x © Ze 2 ) 



We compute the series <£?($, A) = £>(t>i, t> 2 )<it>i<it> 2 where 



e 2«7r(niDi+n 2 f2) 



(2iTTni)(2iTrn 2 )(2iTT(n 1 + n 2 )) 

ni^O, n 2 ^0,ni+ri27^0 

The distribution £>(<3>,A) is piecewise polynomial and periodic with 
respect to A = Zei + Ze 2 . It is thus sufficient to write the formulae of 
B(vi, f 2 ) for < V\ < 1 and < t> 2 < 1 (see Figure HJ). 



„( \ _ f + -2u 2 )(fi - l + i> 2 )(2t>i -f 2 ), if«i<w 2 

-1(^-2^)^ -1+^(2^-1-^), if Vl >va. 

We remark that 3£>(^i,t> 2 ) is the symplectic volume of the moduli 
space of flat SU(3)-connections on a topological torus with one marked 
point v = v\H ai +t> 2 if a2 where H ai and H a2 denote coroots associated 
to simple roots {ai, a 2 } of SU(3). 

Example 2.4. Let V = Re 1 © Me 2 with lattice A = Ze x © Ze 2 . Let 
$ = [ei, e 2 , e x + e 2 , ei — e 2 ]. We write v E V as v = V\t\ + f 2 e 2 . We 
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compute the series £>($, A) = B(vi,v 2 )dvidv 2 where 



g2i7r(ni«i+n2i)2) 



B{Vl,V2) ni^Lez (2inn 1 )(2inn 2 )(2i7r(n 1 + n 2 ))(2i7r(n 1 -n 2 )) 

ni^O, n2y^0, ni+n,2^0, n\— n2^0 

In the region t> i — v 2 < 1, v 2 < 0, v± + v 2 > 0, we get 

1 

B(v 1 ,v 2 ) = -^2(2^1 - 1)K - 1 - v 2 ){vi +v 2 ). 
In the region v± > v 2 ,v 2 > 0,vi + v 2 < 1, we get 

B(ui,u 2 ) = ^v 2 {2v 1 - l)(^i - 1 + v 2 ){v 1 -v 2 ). 

o 

Similar computation in the region v\ > v 2 , v 2 > 0, v± + v 2 < 1 gives 
B(v 1 ,v 2 ) = ^i(2v 2 - l)(^i - 1 + v 2 )(vi -v 2 ). 

3. Recurrence relations 
For an element in $, we associate two lists of vectors as follows: 

• We consider the list $—{0} in V and the corresponding distribution 
B($ - {0},A) on V. 

• Consider the vector space Vq :— V/ < >, let p denote the 
projection V — > Vq. We denote the image under p of the lattice A in Vq 
by A . The dual space U of the vector space V is the hyperplane H^. 
The dual lattice to A is the lattice T = {7 G T | (7, 0) = 0}. Consider 
the list $0 of elements of Vq consisting of the images of the elements 
in $ — {(f)}. Observe that if $ contains with multiplicity greater 
than 1, then <l> contains the zero vector and consequently £>($ , A ) is 
identically zero. 

If D is a distribution on Vo, we denote by p*D the distribution on 
V "constant in the direction 0": if D — b(vo)d\ vo, then we define 
p*D = b{p(y))d^v. Thus p*B(§o, A ) is a distribution on V. We remark 
that, for any G we have the following equality of sets 

(3.0.1) r reg ($ - {0}) = (r rcg ($ - {0}) n {0 = o}) u r rcg ($) 

where the union is disjoint. 

The main remark of this section is the following recurrence relation 
for the distribution £>($, A). 

Proposition 3.1. Let G Then we have 

(3.1.1) c^B(<t>, A) = £($ - {0}, A) - p*B(%, A ). 
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Proof. We fix the measures d^v and oIa v and we identify £>($, A) and 
i3($o, A ) to generalized functions. Differentiating A) in the sense 
of generalized functions, we get 

(2iTTV,~f) 

W) A) ( ,)= y: n ^r—- 

T er reg (*) 1 W-{4>} W>7> 

e (2i7rt;,7} e (2i7rt),7) 

~ 7er re ^- { , }) IW W < W, 7) " 7erreg($ _i; }) , (7 ,, )= o IW W < W, 7) ' 



The last term is constant on the line v + M.<f> and identifies with 
p*B(%,A ). □ 

4. Hyperplane arrangements and generalized series 
We generalize the setting of Bernoulli series. 

Here we assume that the list $ in Vq does not contain the zero vector. 
Then each in $ determines an hyperplane — {u G U : (u, 0) = 0} 
in U. Let 

H($) = {H^ G $} 

be the set of hyperplanes determined by <3>. We denote the closed subset 
U^^Htp of U by the same notation "H($). When $ is fixed, we denote 
"H($) simply by "H, and its complement in U by 

We denote by S(V) the symmetric algebra of V and identify it with 
the ring of polynomial functions on U. We denote by R-u the ring of 
rational functions on U regular on Uu, that is, the ring generated by 
the ring S(V) of polynomial functions on U together with inverses of 
linear forms G $. 

The set r rcg ($) depends only on "H, thus, we shall also denote it by 
r reg ("H). A function g G R-h is well defined at 7 G T reg (H). 

Definition 4.1. If g G Rn, we define the distribution B(H, A, g) on V 
by 

7er rcg («) 

It is easy to see that the above series converges in the space of dis- 
tributions on V. The Bernoulli series £>(A, <£>) is the special case of 
B(n,A,g) with 4 = ^. 

Example 4.2. Let A = Ztu, let U = {0}, and g = 1. Then 

^,A l3 ) = ^e" = -l+5 A 



Multiple Bernoulli series 



13 



Observe that if Ai C A 2 , then B(%, A 2 , g) is obtained from B(H, A\,g) 
by averaging over A 2 /Ai: 

(4.2.1) B(H,A 2} g) = ^ t(X 2 )B(H, A u g). 

A 2 eA 2 /Ai 

Let cf) G $, then we can associate to (f) the following two arrange- 
ments: 

• W = H\ H^. 

• %q = {H n H^, H G "H'}, the trace of the arrangement H' on H^. 

Clearly a function g in restricts to the hyperplane in a ra- 
tional function g lying in Ru - Thus B(T-L , A , go) is a distribution on 

h; = v/r<p. 

We have the following recurrence relation for the distribution £>("H, A, g) 
associated to an element g G Ru 1 ■ 

Proposition 4.3. // g G R-n> , then 

B(H, A, 0) = B(H', A, g) - P *B(U , A , <? ). 

Proof. From the equality (13.0. ip . we see that the elements of r reg ('H') 
that are not in T ieg (7i) can be identified with the elements of r reg ('Ho)- 

□ 

5. PlECEWISE POLYNOMIAL BEHAVIOR 

For completeness we reprove here that the distribution i3($,A) is 
piecewise polynomial when $ generates V. In fact, we prove the piece- 
wise polynomial behavior of the series B(H, A, g) when g belongs to a 
particular subset of Ru which we will shortly describe. 

Suppose $ generates V. A subspace of V generated by a subset of 
elements of $ is called ^-admissible. A ^-admissible hyperplane will 
also be called a wall. Let ^ a ff(^ > ) be the set of ^-admissible hyperplanes 
in V together with their translates with respect to A. An element 
W G H a (i($) will also be called a (affine) wall. An element v G V is 
said to be regular if v is not on any affine wall. We denote by Keg,aff 
the open subset of V consisting of regular elements. We denote by 
T($,A) the set of connected components of V Ieg;SC g. An element r of 
T($,A) is called a tope. By definition, topes only depend on A and 
the arrangement and not on $ itself; thus we will denote the set 

of topes indifferently by T($, A) or T("H, A). 

Suppose Ai C A 2 . Then, topes corresponding to the system ($, A 2 ) 
are obtained by translating topes corresponding to ($, Ai) by elements 
of A 2 and taking their nonempty intersections. 
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f2 




Figure 5. T($, Zei © Ze 2 ) versus T($, Zei © Z^i±M) 



Example 5.1. Let V = Mei©Me 2 , $ = [ei,e 2 ] and A = Zei©Z^±^. 
Let $ = [ei,e2]. Then, the topes in T($,Zei © Ze 2 ) gives a paving 
of V by squares (complement of bold black lines in Figure [5]), and the 
topes in T($, A) are obtained by subdividing the squares into 4 equal 
squares. 

Definition 5.2. A function / on V Tegta s is called piecewise polynomial 
with respect to H a g($) and A if / coincide with a polynomial function 
f T on each tope r in T($, A). 

A distribution D is called piecewise polynomial with respect to %($) 
and A (in short (7i, A), or equivalently ($, A)) if it is given by integra- 
tion on V^eg, a ff by a piecewise polynomial function. The space of piece- 
wise polynomial distributions with respect to (7i, A) is invariant under 
translation by A. More generally, if A x C A 2 , and D is piecewise poly- 
nomial with respect to ("H, Ai), then t(A 2 ).D is piecewise polynomial 
with respect to (K,A 2 ) for any A 2 G A 2 . 

The condition for a distribution b to be piecewise polynomial is 
stronger than the condition that the restriction of b to any tope is a 
polynomial density. For example the S function of the lattice A restricts 
to on any tope r, but is not a piecewise polynomial distribution. 

Let G $, and consider the two arrangements %' and Ho associated 
to 4> as in the previous section. If / is piecewise polynomial for (%', A), 
then / is piecewise polynomial for (H, A). If / is piecewise polynomial 
for (Wo, A ), then p*f is piecewise polynomial for (J-L, A). 

We now prove that B(H,A,g) is piecewise polynomial with respect 
to (H, A) when $ generates V and g is in some subspace of Ry_ that 
we describe now. 

We may assume that all equations = of the hyperplanes in % 
lie in A, we can always achieve this by taking an appropriate multiple 
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of 0. Thus, for what follows, we may assume that elements of $ are in 
fact in A. 

We denote by 23($) the set of subsets of r linearly independent 
elements of $. In other words, an element of *B($) is a basis of V 
extracted from $. 

Suppose L is a sequence of elements of $ (possibly with multiplici- 
ties) generating V. Define 

a function in R^- 

Since 9(L) will change by a scalar multiple when elements of L are 
scaled, we may define the following space which depends only on %. 

Definition 5.3. Let G% be the subspace of Ru generated by all ratio- 
nal functions of the form 9(L). 

We recall the following description of Gy_. 

Lemma 5.4. Any 9(L) may be written as a linear combination of el- 
ements 9(a,n) = Hl 1 nr where a := [a^, . . . ,a ir ] G *B($) is a basis 

extracted from $ and n = [ni, n 2 , . . . , n r ] is a sequence of positive in- 
tegers. 

Proof. By induction on the number of elements of L, we need to prove 
that the assertion holds for rational fractions of the form 9(a,n)-^ r 
with a = [ax, ... , a r ] is a basis of V. We write a = Y^=i c % a %- Using 
the relation, 

at \J_ a 1 V - — 

la ' nJ a N a? 1 ■ ■ ■ a? a N ^ .^> n Q a? 1 • • • a? 4 " 1 • • • a^ 1 ' 

we decrease the number |n| = n\ + • • • + n r . When one of the rij in 
the sum becomes 0, the corresponding term is of the required form 
associated to the basis <jj = o U {a} \ □ 

Proposition 5.5. If g G G%, the distribution 

B(n,A,g)(v)= 9{l)e 2i ^ v) 
7er rcg (w) 

is a piecewise polynomial distribution (with respect to the system (TL, A) ). 

Proof. We prove the proposition by induction on the number of ele- 
ments in %. Using Lemma I5.4[ it suffices to prove the proposition for 
g of the form g = 9(a,n) for cr G *B($) and n a sequence of positive 
integers. 
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We first assume that $ consists of independent elements cti, 0*2, ■ ■ ■ ,a r 
possibly with multiplicities. Consider an element g = 9(L) of Gy for 

L = [oci, . . . , ocx, a 2 , ■ ■ ■ , a 2 , • • • , a r , . . . , a r ] 

where a\ appears with multiplicity ki in $. 

Let A' = Zai © Za 2 © • • • © Za r ; clearly A' is a sublattice of A. We 
choose coordinates t = J2i^i a i- Then B(T-L,A',g) = B{t) with 

(5.5.1) B(t) = (— l) r TT -^-rB(ki, U - [U]). 

The function -B(t) is a polynomial function on each parallelogram trans- 
lated from the parallelogram 5^[ =1 [0, l]aii. By equation (17.5.ip . the 
distribution B(T-L,A,g) is obtained by averaging B(T-L,A',g) over A/A'. 
Thus B(%,A,g) is piecewise polynomial with respect to (%, A). 

We now consider the general case, where the cardinality of % is 
greater than r. In this case there exists G $ with the property that, 
for the hyperplane arrangements %' and Hq associated to 0, we have 
g G Gyi and g (the restriction of g to (f) = 0) is in G^ . Then, using 
Proposition I4.3[ which states 

A, ff ) = B(H', A, (/) - p*B(Uo, A , ff0 ), 

we conclude by induction that B(H,A,g) is piecewise polynomial with 
respect to (H,A). □ 

The Bernoulli series B($,A) is equal to BCH,A,g) with g = n 1 • 

it is an element of for we assumed that $ spans V. Thus, by 
Proposition I5.5[ we immediately obtain that B($>,A) is a piecewise 
polynomial density. 

Corollary 5.6. For any / G Ry, the distribution B(H,A,f)(v) re- 
stricts to a tope polynomial density. 

Proof. Let / G R%. We can write / as Pg where P is a polynomial and 
g G Gy- Then the distribution B(%, A, f)(v) is obtained by applying 
the differential operator P(d v ) to the distribution B(H, A, g){v). 

This differentiation is in the distribution sense so that it may produce 
distributions supported on admissible hyperplanes, but on an open tope 
r, we obtain a polynomial density. □ 

Definition 5.7. Given a tope r in T($, A), we denote by Ber($, A, r) 
the polynomial function on V such that the restriction of B(<&,A) to r 
coincides with the restriction of Ber($, A, T)(v)dxv on r. 



Multiple Bernoulli series 



17 



By the above proof, we see that the polynomial Ber($,A,r) is of 
degree equal to the number of elements in $. 

The fact that £>($, A)(t>) is a periodic distribution on V implies im- 
mediately the following periodicity formula. For any A G A and v G V, 

(5.7.1) Ber($, A, r + X)(v + A) = Ber($, A, r){v). 

If v is a connected subset of V contained in the open set of regu- 
lar elements, we denote by Ber($, A, v) the polynomial Ber($, A, r(z/)) 
where t{v) is the unique tope containing v. 

Let r G T($, A) and <p G $. If u Q G V = V/ < > is the projection 
of v G r, then w is not on any affine wall in Vq. Indeed the reciproc 
image of an affine wall in Vq is an affine wall in V. We denote by Tq 
the unique tope in Vq containing the projection of r. 

Equation (13.1. ip implies the following relations. 

If $ — {0} generates V, then 

(5.7.2) Ber($, A, r) = Ber($ - {0}, A, r) - Ber($ , A , r ). 
If $ generates V, but $ — {</>} does not generate V, then 

(5.7.3) fyBer($,A,r) = - Ber($ , A , r ). 

Remark 5.8. By using reduction to independent variables and the 
explicit formula (15.5.11) . we obtain also a way to compute Ber($, A, r). 
This can be applied not too painfully when the number of elements in $ 
is small. However, the residue formula due to A. Szenes [12] to compute 
Ber($, A, r) is very efficient when $ is large, provided the dimension of 
V is relatively small. We will give examples of computations of volumes 
of moduli spaces using Szenes formula in a next article. 

6. An Euler-MacLaurin formula 

This section is independent of the rest of the article. 

Assume that $ generates V. Using the Lebesgue measure associ- 
ated to the lattice A, we identify i3($, A)(t>) to a piecewise polynomial 
function on V. 

Let us denote by 71 the set of ^-admissible subspaces of V. Then 
s = V and s = {0} are the maximum and minimum elements of the 
partially ordered set 71. If s is a ^-admissible subspace of V, we denote 
by $ \ 5 the sequence of elements in $ not lying in the space s. 

The projection of the list $ \ s on V/s will be denoted by $/s. 
The image of the lattice A in V/s is a lattice in V/s. If $ generates 
V, $/s generates V/s. Using the projection V —> V/s, we identify 
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the piecewise polynomial function B(<&/s, A/s) on V/s to a piecewise 
polynomial function on V constant along the affine spaces v + s. Then, 

r rcg ($/s) :=rn[/ rcg ($/s) 

is the set of elements 7 G T satisfying (7, s) = for all s G s and 
(7, 0) 7^ for all G $ \ s. 

We lift functions on V/s to functions on V by the canonical pro- 
jection. Thus £>(<£>/s, A/s) is the function on V given by the series 
(convergent in the sense of generalized functions) 

„^ e 2tir(«,7) 

This function is periodic with respect to the lattice A, piecewise poly- 
nomial on V (relative to ($,A)) and constant along v + s. We denote 
it simply by £>($/s) leaving its dependence on the lattice A implicit. 
If s = V, the function £>($/s) is identically equal to 1; if S = {0}, then 
we obtain the multiple Bernoulli series £>($,A). 

Theorem 6.1. Let f be a smooth function on V , rapidly decreasing 
with rapidly decreasing derivatives. Then, 

E/( A ) = E(- 1 ) l * Vs| / B w*)( v )( n woo**- 

AeA sen v <^e<i>\s 

Remark 6.2. The term corresponding to s = V in the above sum 
gives the term j v f(v)dv. Thus we may also write the formula as 

£/(A)- / f(v)dv = J2(-^ sl [ B($/ S )(v)(H d+)f{v)dv. 
AeA Jv b^v Jv 0e*\s 

All the sets $ \ s entering in this formula are 'long', that is, their 
complement in <3> do not generate V. In particular they contain a 
'cocircuit'. This formula has been used in [16J to obtain a formula for 
the semi-discrete convolution with the Box Spline. 

Proof. Let 

f(y) = [ e 2 ^f(x)dx. 
Jv 

By Poisson formula 

AeA 7 er 

We group together the terms in T belonging to s -1 for s G TZ. More 
precisely, the lattice T is a disjoint union over the s G TZ of the sets 

r reg ($/s) = {7 G S 1 n T\(<f>, 7) ^ for all G $ \ s}. 
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Now in the generalized sense 

E fa) = I E ^^f(v)dv 

7 Gr reg (*/«) Jv ^v^/s) 

d^ s B(^/s)(v)f(v) 

'v 

and we obtain the statement in the theorem. □ 

7. Wall crossing 

In this section we again assume that $ generates V. Under this 
assumption, we compare the polynomials Ber($,A,r) associated to 
two adjacent topes of T($,A) separated by an hyperplane W. We 
remark that due to the periodicity property of i3($, A) it suffices to 
consider jumps over an hyperplane W passing through the origin. 

If D 1 and D 2 are two distributions on V with supports Si and S 2 
with the property that for any v G V the intersection of v — Si and 5*2 
lies in a compact set, then the convolution Di * D 2 is well defined. 

We recall the definition of multispline. Let X = [vi,v 2 , ■ ■ ■ ,v m ] be 
a sequence of non-zero vectors in V. We will first consider the case 
where X spans a pointed cone. The multivariate spline T(X) is the 
tempered distribution defined on test functions / by: 

POO POD m 

(7.0.1) (T(X)\f)= ••■/ f(52t i Vi)dti---dt m . 

Jo Jo i=1 

If X spans V, we may interpret T(X) as a function on V supported in 
the cone C(X) generated by X. This function is piecewise polynomial. 
If v G X, then d v T(X) = T(X - {v}). When X is the empty set, 
T(X) = 5 . 

We now consider the case where the elements of X do not necessarily 
lie in a half-space. We introduce a polarization of X given by a vector 
u in U. Let u G U be a vector that is nonzero on all elements of X. 
We will then say that the vector u is polarizing for X. Divide the list 
X into two lists A and B, the lists of positive and negative vectors on 
u respectively. We then define 

T{X,u) = {-l)WT([A,-B\). 

Example 7.1. With the notation of Example U.2\ 

, ti^t if t > 0, , „ { if t > 0, 



(fc-l)! 
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We return to our set up. Let $ be a sequence of nonzero vectors in 
V, spanning V. Let W be a ^-admissible hyperplane. Let E G T be 
an equation of this hyperplane, where E is a primitive vector in T; this 
fixes E up to sign. The lattice A is fixed, and we write simply dv instead 
of d\v . Similarly we denote by dh the density determined by Anil 7 . As 
E does not vanish on any element of $ \ W, we may define T($ \ W, E) 
as above; it is a distribution supported on E > 0. Let p be a polynomial 
density on W. Then, the convolution p*T(<&\W, E) is well defined and 
it is supported on E > 0. Similarly, p * T($ \ W, —E) is supported on 
E < 0. It is easily proven (see pQ) that p*T(<S>\W, E)-p*T(<S>\W, —E) 
is given by integration against a polynomial density. We thus define the 
polynomial Pol(p, $\W,E) by the equation 

Pol(p, $ \ W, = p * T($ \ E) - p * T($ \ W, -E). 

The following properties of Pol(p, <&\W,E) follow directly from the 
above equation. 

Lemma 7.2. Let * = $ \ 

^ajLei ^ G T/ien, 



(^cj // l^l > 1, then the restriction of Pol(p, E) to W vanishes of 
order — 1. 

The following one dimensional residue formula for Pol(p, <&\W,E) 
is given in pQ. It is useful in computing the convolutions. We write 
p(h) = f(h)dh where / is a polynomial function on the hyperplane W. 

Lemma 7.3. Let P be a polynomial function on V extending f . Then, 
for v G V, 




Pol(p,$\W,E)(v) = Res 



2 = 




P(d x ) ■ 



U^\w(^^ x + zE ) 



) 



x=0 



) 



Theorem 7.4. Let T\ and T2 be two adjacent topes in T(<&, A) separated 
by the hyperplane W defined by E with (v,E) > for any » 6 tj. 
Denote by r 12 the tope in T($ fl W, A fl W) containing rj" fl in its 
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closure. Let Ber T12 := Ber($ n W, A fl W,Ti 2 )dh be the polynomial 

density on W determined by T\ 2 . Then, 

(7.4.1) 

(Ber($, A, n)-Ber($, A, r 2 ))dv = Ber T12 *T($\W, E)-Ber T12 *T($\W, - 

Remark 7.5. Formula (17.4. ip is very similar to jump formulae for vol- 
ume of reduced spaces in Hamiltonian geometry. Indeed if fi : M — > t* 
is a proper moment map associated to an Hamiltonian action of a torus 
T, then the set of regular values of fi is the complement of a certain 
number of affine hyperplanes. On each connected component, vol- 
umes of reduced spaces M red {v) := /i _1 (t>)/T are given by polynomial 
functions of v. When crossing a wall, the variation of these polyno- 
mials follow the same jump scheme as in equation (17.4=. ip : they are 
determined by a polynomial volume function associated to a smaller 
Hamiltonian manifold Mo and weights of the normal bundle of Mo in 
M ([§]). 111 particular, when the sequence $ is comprised of positive 
coroots of a compact connected Lie group G with multiplicity 2g — 1 
and A is the coroot lattice of G, the polynomials Ber($, A, r) describe 
(up to some normalization) the symplectic volume of the moduli space 
of flat G-connections on Riemann surface of genus g with one bound- 
ary component, around which the holonomy is determined by v. These 
moduli spaces are reduced spaces M re( i(v) of an Hamiltonian action. 

Proof. We will first verify the claim for the case where there is only one 
vector in $ that is not contained in W . 

Let A = A n W. We consider the lattice A b = A © Z0. By formula 
fl2XID, 

(7.5.1) B($,A)(v)= J2 B(*,M(v + *i)- 

XjeA/At 

For t in a small neighborhood of zero, let t\ (respectively r%) denote 
the tope containing the open set of v — h + tcf) for h G W lying in a 
relatively compact open subset of t\i and t > (respectively t < 0). 

We may express a representative of a non-zero Xj G A/A& as Xj = 
hj + tj(f> for hj G A and tj Z. As the lattice A b is product of lattices, 
we have 

A 6 ) (h + t<p) = B{$ n W, A ) (h)(-t + [t] + -)dt. 

Observe that the jump in the function £>($, A b )(v + Xj) = B($, A b )(h + 
hj + (t + tj)4>) as t changes sign in a small neighborhood of zero is 
precisely zero for the nontrivial representative since tj is not integral. 
Thus the only contribution to the jump comes from the trivial Xj in 
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the sum of equation (17.5.11) . We get 
Ber($, A,Ti)(v) - Ber($, A, t 2 )(v) 



' Ber($nW,A )(/i)(H+ I)). 



(E,cf>) " ' 2 ^ 2' 

The convolution product in this case is just the product in coordi- 
nates so that 

Ber($ n W } A , r 12 ) * T({0}, E) - Ber($ n W, A , r 12 )(/i) * T({0}, -£) 

is equal at the point (h, t) to Ber($nW / , A , Ti 2 )(h) and hence we obtain 
the claimed formula. 

Now consider the case where there are several elements of $ that do 
not lie in W. Let be a vector in <3> \ W. Let $' = <& — {0}; $' still 
generates V and $'D W = ^HW^. Equation (13.1.11) implies that in this 
case A) is continuous on W: indeed the derivative in the direction 
is a piecewise polynomial function. 

Let t[ and r 2 be the topes of $' containing n and r 2 respectively. 
They are adjacent with respect to and Ber T12 = Ber Tl2 . Using equa- 
tion (15.7. 2p . we have 

d+ Ber($, A, t x ) - Ber($, A, r 2 ) = Ber($\ A, r[) - Ber($', A, r' 2 ). 

Indeed the topes T\ and r 2 give the same tope r under projection onto 
V = V/ < >. 

By Lemma [7.21 part (a), 

^Pol(Ber T12 , $ \W,E) = Pol(Ber r «, $' \ W, £). 

Denote by Leq($) the left hand side, and by Req($) the right hand side 
of equation (17.4. ip . By induction, we have 9^(Leq($) — Req($)) = 0. 
Thus, the polynomial function is constant in the direction of K.0. The 
left hand side vanishes on W by the continuity of -B(<3>, A) on W. Hence 
the claim. □ 

We now demonstrate the theorem with various examples. 

Example 7.6. Recall the data of Example 12.21 Let T\ and r 2 be two 
adjacent topes defined by inequalities < t < 1 and — 1 < t < 
respectively. By Theorem 17.41 and Example I7.1[ 

a-k-l 

Ber($, A, n)(tu) - Ber($, A, r 2 )(tu) 



(k-l)V 

which is indeed equal to —hB(k,t) + hB(k,t + 1) as it can be seen 
from the explicit expression of A)(tu) in Example 12.21 



Multiple Bernoulli series 



23 




Example 7.7. Recall the data of Example 12.31 Let T\ and r 2 be the 

two adjacent topes separated by the hyperplane W = M(ei + e 2 ) (see 
figure HD^a)). Then E = — e 1 + e 2 . 

We express v = v\&\ + t> 2 e 2 G V as v = vi(ei + e 2 ) + (u 2 — Vi)e2 
and x G U as x = Xie 1 + x 2 (— e 1 + e 2 ). Using Example 12. 2^ with 
$ n W = ei + e 2 and Anlf = Z(ei + e 2 ) at ^i(ei + e 2 ), we have 

Ber T12 = Ber($ n W, A n W, 7i 2 )(vi(ei + e 2 )) 
= 5 - «i 

In the above coordinates of U, the operator Ber T12 (9 x ) = | — (9 Xl under 
the identification P(d x )e ( - a ' x ^ = P(a)e ( ~ a ' x K Then, 

Pol(Ber- d>, £)(,) = Res z .o ((bct-^) • n^g^ 



x=0 

Res z=0 ( ( (5 - • {X1 _ X2 _ Z)(X2+Z) 

e vix 1 + (v2-v 1 )z 



Res z=0 m| -dx x 



x=0 



- |(1 - -u 2 )(«i -f 2 ), 

which is indeed the jump Ber($, A, T\){v) — Ber($, A, r 2 ) (v) as it can 
be seen from the explicit expression in Example 12.31 

Example 7.8. Recall the data of Example 12.41 

(a) Jump over the wall W = WLe±: Then E = e 2 , <3> HW = {ei} and 
A n W = Zei (see figure EJo)). 



Poller- $,£)(«) = Res 2 = (W-(^) • 

= res, =0 ((-^ + |) • - x] ^- + -— 



x=0 . 



2+z)(x 1 +X2+z)(xi-X2-z) J _ Q 

H( 2u i - 1) 



2 = 

1 
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which is indeed the jump Ber($, A, tq)(v) — Ber($, A, Ti){v) as it can 
be seen from the explicit expression in Example 12.41 

(b) Jump over the wall W = R(ei + e 2 ): Then E = — e 1 + e 2 . 
We express v = V\e% + t> 2 e 2 G V as v = v\(ei + e 2 ) + (f 2 — i>i)e 2 
and x e U as x = X\e l + X2 (— e 1 + e 2 ). Using Example 12. 2[ with 
$ n PU = ei + e 2 and A D PU = Z(ei + e 2 ) at v 1 (e i + e 2 ), we have 
Ber 7 " 23 = Ber($ nW,AnW, T 23 )(v 1 (e 1 + e 2 )) = § - v v Then, 

Pol(Ber™, B)(„) = Res„ (((| - 9x0 ■ ^^l-Sl, ) 

= -!(«! - 1 + t> 2 )(t>i -t> 2 ) 2 , 

which is indeed the jump Ber($, A, t 3 )(v) — Ber($, A, t 2 )(v) as it can 
be seen from the explicit expression in Example 12.41 



8. The affine case 

This section generalizes previous results to the affine case. Results 
proven here are not needed for the following section. 

Let $ = . . . , 0tv] be a list of elements of Vq and let z = [zi, z 2 , ■ ■ ■ , z N 
be a list of complex numbers. We consider the augmented list $ := 
[[0i, Zi],..., [</>n, zn]\ and define 

r reg ($) = rrw rcg (<J) 

where 

£/re g ($) = {ueU\ (0J, u) + Zj ^ for all j}. 
Definition 8.1. The affine multiple Bernoulli series is the distribution 

S($ ' A)= T, n ^v 9 . 777 r- 7 ^- 



The distribution A) has the following properties, similar to its 
nonafline counterpart: 

• If $ is the empty set, then £>($,A) is the 5-distribution of the 
lattice A. 

• If Ai C A 2 . Then i3(<3>, A 2 ) is obtained from B(&, Aj) by averaging 
over A 2 /Ai: 

(8.1.1) iB($,A 2 )= ^ t(A 2 )i3(<l,Ai). 

A 2 GA 2 /Ai 
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In the special case Zj = (4>j,z) for z G Uc, it is more natural to 
consider the distribution 

g2i7r(?;,7+z) 

Eto(*, A, *)(*)= nN - 77 - r - 

7er;(0J , 7+2 )^o IWM^T + z) 

Clearly, 

Eis($,A,z)(v) = e 2i7T{v ' z) B($, A)(v) 

for $ = 2)], • • • \4>n -i ((j>N, z})]. If z is regular, that is (cj),j}+z ^ 

for all G $, then Eis($, A, z)(t>) defines a distribution of v with 
coefficients meromorphic functions on Tq = Uc/T which is studied in 

®- 

Example 8.2. Let A = Zu, and let = [[u, z], [cj, z], ... , [w, z]] where 
[a;, z] is repeated k times. If z is integral, we simply have 

B(§ k ,A)(tu) = e~ 2i7rzt B(§ k ,A)(tu). 

If z is not integral and k = 1, then, using Lemma 16 of [3], 

^liimt g([t]—t)2iirz 



(8.2.1) B(* u A) (to;) = J2 9 . , ■ A = T 

li'Kva + z) 1 — e 



which is an analytic function of t in each tope. 

If k > 1, z not integral, and < t < 1, we use the residue theorem 
for the integral 



e~ tu 

— -du 



I\u\=r (2iTTZ-u) k (l-e—) 
which tends to when R tends to infinity. Then 

B(i>fe, A) (tu) = - ReS u=2 inz (2inz-u)*(l-e-v) du 
-2inzt r> j 



Thus, we see that $($/-, A)(to) is a product of an exponential function 
of t and a polynomial in i. In particular, in the interval < t < 1, it is 
an analytic function of t. For example, for k = 2 and < t < 1, we get 



-H-KZt 



g($ 2 ,A)M= „_»■ (* + 



^ g — 2i7T2: \ g2«7TZ J 

8.1. Recurrence relations. In the affine case the recurrence relation 
( 13.1. ip is slightly modified. 

Let = [0, z] be an element of $. We consider two cases. 

• Suppose there exists ^ z G T such that 

(8.2.2) (0, 7z ) + z = O. 
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Then, we may express 7 G r reg ($ — {0}) satisfying (7, 0) + z = as 
7 = 7' + 7 2 - Clearly, (7', 0) = 0. 

We consider the system 

$ = [[0 i; Zj + (0^7,)], <t>j e $ - {0}] 

in Vo = V7 < >. The sum 

e (2i7ro,7 / ) 



E 



is constant in the direction of and identifies with i3($ , A)(i>). Hence, 

we get the following recurrence relation. 

(8.2.3) 

(8+ + 2iTrz)B($,A)(v) = 8($ - {0}, A)(v) - e 2 ™^B(® ,A )(v). 

• If there does not exist 7^ satisfying the relation (I8.2.2p . then 
r r eg(^ — {0}) = r rcg ($), and the equation f!8.2.3|) reduces to 

(fy + 2iirz)B@,A)(v) = £?(<£ - tf},A)(i)). 

8.2. Piecewise exponential polynomial behavior. We consider 
= [0, 2] G $ with 7^ 0. Consider the complex hyperplane Hi : = 
{u G Uq '■ (u, 0) + z = 0}. Consider the set 

H = U$) = {Hi, 0G$} 

of hyperplanes in Uq- We denote by TZ^ the ring of rational functions 
on Uc with poles along %. That is, if S(Vc) denotes the symmetric 
algebra of Vc, identified with the ring of polynomial functions on Uc, 
then IZ-fi is the ring S(Vc) of polynomial functions on Uq together with 
inverses of forms (0, •) + z for [0, z] G $. 

For (7 G 7^ we define the distribution B(H,A,g) on V by 

B(fi,A,g)= 9(l)e 2l7T{v ^d A v, 
7er rcg («) 

where r reg (?^) = r rcg ($), as regularity does not depend on the multi- 
plicity of an element in $. 

We fix G $, and define %' := % \ Hi. For g G 72.^,, we compare 
B(%A,g) and B(W,A,g). 

Similar to the nonaffine case, for a fixed G we define T-L' := H\H^ 
and T-^o to be the collection of affine hyperplanes H fl Hi for those 
H G "H not parallel to Hi, that is, for H associated to [<f>j, Zj] G $ 
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with (f)j 7^ <p. The collection "Ho is a collection of affine hyperplanes in 
the affine space Hi. 

We consider two cases: 

• There exists 7 Z G T lying in Hi. Thus ((ft, j z ) + z = 0. Let H be 
the real hyperplane with equation = 0. If K G Hq, then if — {7^} is a 
complex hyperplane in (Hq)c- Let "Hg be the collection of hyperplanes 
K — {7 2 } with if G "Ho- Then, for g G i?^/, we define go(u) := <3 l (tt + 7.z) 
lying in TZ^. Let Vo = V/M(p, and A the image of A in V . 

It immediately follows from the set theoretic partition in the proof 
of Proposition 14.31 that: 

Lemma 8.3. If g G Rfir, then 

B{U, A, g) = BCH', A, g) - e 2 ^ v ^p*B(n z , A , g ). 

• In the case that there does not exist any 7 2 G T lying in Hg and 
satisfying Equation (18.2.21) . we have 

B(H,A,g)=B(H',A,g). 

For a fixed $ we will denote the list of vectors coming from the 
first component of the pairs in $ by $. Suppose that the vectors in $ 
associated to $ span V. Let Gfi denote the subspace of TZ^ generated 
by functions of the form 

0(L)(x) = 1 

UaeL {U 1 x)+Z a 

where L is a list of vectors coming from $ generating V. 

We call a function that is a sum of products of exponential functions 
and polynomial functions an exponential polynomial. 

We will say that a locally L 1 function / is piecewise exponential poly- 
nomial, if there exists a decomposition of V in a union of polyhedral 
pieces Cj such that the restriction of / to Cj is given by a exponen- 
tial polynomial formula. We then say that the distribution f(v)d/^v is 
piecewise exponential polynomial. 

Proposition 8.4. If g G Q-n, then B (ft., A, g) is a piecewise exponential 
polynomial distribution. 

Proof. We use the same line of argument as in the proof Proposition 
15. 51 As before, we scale the denominator of g = 9(L) such that all 
a G L lie in the lattice A. In the case that L has independent elements, 
B(l-L, A, g) can be written as a product of exponential polynomial func- 
tions £>(<3>fc,A), whose expression changes whether the (scaled) z are 
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integral or not. The expression for both cases is given explicitly in 
example 18.21 and they are piecewise exponential polynomials. We then 
use the averaging formula (18.1. ip . 

In order to reduce the general case to the case of independent vectors 
we use an analogue of Lemma I5.4[ and in the case that same a with 
distinct z appears in $, we use the relation 

1 1111 

+ 



(a + Z\){ot + z 2 ) Z\ - z 2 (a + Zi) z 2 — Z\ (a + z 2 ) ' 

We then get the claimed property of B{7i, A, g) by induction using 
Lemma 18.31 □ 



The above proposition for $ = [[(pi, z\], . . . , [(f) n, zn]] and 

g{x) 



nf =1 2z7T((0 J ,x) + 



Lj=l aun i z jJ 

gives: 

Corollary 8.5. If $ associated to $ generates V, then B(<&,A)(v) is 
an exponential polynomial function of v on a tope of T($, A). 

Remark 8.6. Using the same proof as above, we see that Eis($, A, z)(v) = 
e 2in ( v,z 'B(&, A)(i;) is an exponential polynomial function of v on each 
tope r in T($,A). Furthermore, when z is regular, the recurrence 
relation simplifies to 

(8.6.1) Eis($, A, z) = Eis($ \ {0}, A, z). 

The system of relations in (I8.6.ip are the relations of Dahmen-Miccelli 
[3]. In particular, on each tope r, we obtain that Eis($, A, z)(v) = 

Ki(v)Fi(z) where Ki(v) are Dahmen-Micchelli polynomials and Fi(z) 
meromorphic functions of z. 

8.3. Wall crossing. Given a tope r in T( ( I ) , A), we denote by Ber(<3>, A, r) 
the polynomial exponential function on V such that the restriction of 
A) to r coincides with the restriction of Ber($, A, r){y)dv on r. 
Let (H(a,z)\f) = J t>Q f(ta) e - 2Mz dt. 

Given a wall W, assume that we have renumber $ so that $ = 
[[0i, zj, . . . , [(p p , z p ), [(f) p+1 , z p+1 ), . . . , [(f) p+q , z p+q \) where the first p ele- 
ments (pk belongs to W and the last q elements 4> p+ j do not belong to 
W. Then, we define the lists 

<JrW := [[(f) 1 ,z 1 },...,[(f) p ,z p }} 

and 

$ \ W '.= [[(f)p+U Zp+l], • • • , [^p+g, Zp+q]]- 
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Let E be an equation for the wall W. We define 

T@\W,E):= H -Hi-fa, - Zi )* J] Hfazi). 

(<f>i,E)<0 (<k,E)>0 

We remark that due to the periodicity property of £>(<!>, A) it suffices 
to consider jumps over an hyperplane W passing through the origin. 
We have, similar to Theorem I7.4[ 

Theorem 8.7. Letr\ andr 2 be two adjacent topes o/T(<&, A) separated 
by the hyperplane W , with equation E. Assume that (v,E) > for 
v G T\. Denote by r 12 the tope in T($ PI W, A n W) containing r{ fl 
in its closure. Let Ber($ fl W, A fl W, T^dh be the analytic density on 
W determined by ri 2 . Then, 

(Ber($, A, t x ) - Ber($, A, r 2 ))dv 

= Ber(<l n W, r 12 ) * T($ \ W,E) — Ber(l> n W, r 12 ) * T(<f> \ W, —E). 

Proof. The proof follows the same line of argument as in the proof of 
Theorem I7.4I For the first inductive step, we are reduced by the same 
argument as in Theorem I7.4I to a product situation of W with the line 
M.(f>. Then we compute explicitly using Formula 08. 2. II) . □ 

Example 8.8. Recall the data of Example I8.2I Let T\ and r 2 be two 
adjacent topes defined by inequalities < t < 1 and — 1 < t < 
respectively. By Theorem 18. 7\ 

Ber(i» 1 , A, Tl )(t) - Ber(i> 1 , A, r 2 ){t) = e~ 2 ' mzt , 

which is also seen from the explicit expression of B(<&i, A)(tw) in Ex- 
ample E3 

9. A DECOMPOSITION FORMULA 

Let A and $ be as before. We do not necessarily assume that $ 
generates V. 

In this section, we express £>($, A) as a sum of distributions A(Q, A, a, (3) 
associated to affine admissible subspaces a and a generic vector (3 in 
V. 

Let us start the construction of the distribution A(&, A, a, 0). 

Let s be a ^-admissible subspace of V. Then $ fl s generate s, and 
A fl s is a lattice in s. Let r be a tope in T($ fl s, A fl s). We can 
then consider the distribution £>(r)(s) := Ber($ Pis, A fl s,r)(s)ds. It 
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is a polynomial density on s. We still denote by B{j) this distribution 
considered clS du distribution on V: 



Let A G A. Then a := A + s is an affine ^-admissible subspace of V. 
We say that a is of direction 5. By definition, a tope r of o is such that 
r — A is a tope in s. We define £>($ Hs, r) as a distribution supported 
on a by the formula 



We remark that the definition of 5($ R s, r) above depends only on 
t and not on the choice of A. Indeed, for another A' G A such that 
a = A' + s, X' is necessarily of the form A' = A + A for some A G A Pis. 
Then, 



test(s + A + Ao)S(r-A-A )(s) = / test(s + X)B(r-X- A )(s-A ). 



Using relation ( 15.7.11) . we have B{r — A — An)(s — Ao) = B{r — A)(s), 
hence the independence of the expression. 

For a ^-admissible subspace s, consider an element u G U vanishing 
on s and polarizing for <3> \ s. Then, the multispline distribution T($ \ 
S,u) is well defined. 

Definition 9.1. Let a be a ^-admissible affine subspace of V of direc- 
tion s. Let t be a tope in a, and let u G U be a vector vanishing on s 
and polarizing for $ \ s. Then, we define 



The distribution *4.(<3>, A, o, r, u) is supported on a + m> . It is poly- 
nomial in the direction s. 

Remark 9.2. Choose a direct sum decomposition V = s © r and 
express d e K as t) = s + r for s G s and r G r. If $ is equal to 
$ns©$nt, then the function A(§, A, a, r, w) is, in product coordinates 
(s, r), the product of 5($ Ds, r)(s)(is with T($ \s, w)(r). In general it 
is still possible to express A(&, A, a, r, w)(s,r) as a linear combination 
of product of multispline functions on c and polynomials on s. 

Our main theorem is that £>($,A) can be decomposed as a sum of 
distributions ^4($, A, a, r, u) over all ^-admissible affine subspaces a 
for conveniently chosen r and u. Thus we think of the distributions 
A(§, A, a, r, w) as the basic building blocks of the theory. 





A(<$>, A, a, r, u) := £($ n s, t) * T($ \ 5, u). 
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Choose a scalar product (, ) on V. If W is a subspace of V, or a 
quotient space of V, then W inherits a scalar product. 

Let (3 G V, and let a be a $ admissible affine subspace of direction 
5. We can then write (3 = (3 — f3\ where (3 E a and j3\ 6 s 1 . The point 
Po is the orthogonal projection of (3 on a. Assume (3 generic so that 

• the point (3 lies in a tope t(/3 ) of a. 

• the element f3\ is polarizing for $ \s: (<p, Pi) ^ for all G $ and 
not in s. 

We can then define 



Here the sum is over all admissible affine subspaces a. 

The sum above is infinite. But remark that, given a vector v G 
V by the definition of A(§, A, a, f3), there exists only finitely many 
^-admissible affine spaces a such that A(§,A, a, (3) gives a non zero 
contribution at the element v G V, therefore the above sum is well 
defined. 

For example, if s = 0, then the affine spaces a of direction s are 
reduced to the points A in A and 



We see that A(§, A, {A}, (3) is supported in an affine space A + £ with 
(£, A — f3) > 0. Thus the points v in the support satisfy ||t>|| 2 > ||A|| 2 — 
\\f3\\ 2 . In particular the sum of the distributions 



is well defined. Similar estimates hold for any admissible subspace s, 
when considering the sum over all affine spaces of direction s . 

Remark 9.4. If $ generates V, then V is admissible, and the term 
corresponding to V is the polynomial density Ber($, A, t((3)), with t((3) 
the tope containing (3. The other distributions A(§, A, a, (3) with a ^ V 
are piecewise polynomial densities with support not intersecting r. 

Theorem 19.31 has the following meaning: although the distribution 
£>($, A) is very complicated, it is however obtained by superposing 
simpler functions which are products of polynomials and multisplines. 

Before giving the proof of this theorem we demonstrate the decom- 
position in various examples and state a recurrence relation. 




a 



A^,A,{X},(3) = 6 x *T^,X-(3). 



^^($,A,{A},/3) 
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(a) (b) 
Figure 7. The decomposition of <B($ 2 , A)(tw). 

Example 9.5. Let $ = 0. Then, by definition, A)(v) = E 7 eA* e 2 ™^> v) . 
We would like to decompose this sum as in equation (I9.3.ip . Observe 
that in this case a consists of points of A and any /3 in V is generic. 
Using T(0,/3 — A) = <5o, the decomposition in Theorem 19.31 gives 

7eA* AeA 
which is the Poisson formula. 

Example 9.6. (one dimensional case) Let A = Zu;, and let = 
[u, uj, . . . , u], where u is repeated k times. Then s = 0oys = V = Rcu, 
correspondingly a are reduced to points {A} in A or a = V. Choose 
any (3 = ru G V with < r < 1, it is generic. The polynomial 
Ber($, A, r(/3))(tu) which coincide with £>($ fc ,A)(i) on < t < 1 is 
— jjB(k,t), where B(k,t) is the Bernoulli polynomial. Then, with the 
notation of Example 17.11 Theorem 19.31 gives, 

B($ h ,A){tu)= -±B{k,t)dt+ 8ru,*T{<S> k ,u*) 

In Figure[7]we depict the decomposition of B($2, A)(tu). In part (a) 
we draw the graph of the periodic polynomial — |-B(2, t — [t]), the red 
graph in part (b) is the graph of the polynomial — \B{2 1 1) and lines in 
black correspond to contribution of splines. 

Let us now study the recurrence relations that the distributions 
*4.($,A, a, 0) satisfy. It will be convenient to define A(&,A, a, (3) for 
any affine subspace a, by declaring it to be equal to zero if a is not 
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admissible. If a = A +s where G s, we denote by a/ < > the image 
of the rational space a in V = V/ < >. 

Lemma 9.7. Let G $ and /3 G V. VFe sfe'W denote by (3 the projection 
of (3 onV/<(f)>. 

(i) If <p ^ s, then 

r3^($,A,a,/3)=^($\{0},A,a,/3). 

(ii) If G s, then 

rVl($, A, a, 0) = -4($\{0}, A, a, p)-A{$/ < >, A/ < >, a/ < >, /?). 
In (ii), „4($ \ {0}, A, a, /3) is zero when $ Pis \ {0} does not generate 

3. 

Proof. Part (i) follows from the relation 

^T($\ 5 ,/3 1 )=T(($\{0})\ 5 ,/3 1 ). 

For part (ii) we use equation (13.1. ip on r, which gives 

dtp Ber($Hs, AOs, r) = Ber(($Hs)\{0}, AOs, r)-Ber($ns/ < >, AOs/ < >, r). 

Now suppose g is a polynomial function on s constant in the direction 
of 0. Let X := [t>i, t> 2 , . . . , Ujy] be a sequence of nonzero vectors in V\s 
generating a pointed cone. We denote the projection of Uj to V/ < > 
by V{. Then, 



/»oo /»oo 

Jo Jo 

poo >oo 

= / ' • • / q(v- ^ Uvijdtx ■■■dt N 
Jo Jo 

= q*T(X/ < >){v). 

Putting g = Ber($ ns/ < >, A n s/ < >, r) and X = $ \s, we get 
part (ii). □ 

Proof. We now prove Theorem 19.31 by induction on the number of el- 
ements in $. We assume that the theorem is true for any sublist of 
$. Denote by Req($) the right hand side and by Leq($) the left hand 
side of equation (19.3. ip . 

Let G $. Let $' = $ \ {0}. We recall equation (I3.1.ip . 

fyB($, A) = B{&, A) - < >, A/ < >). 

Let 7?. a fj( < I ) ) be the collection of all ^-admissible affine subspaces of 
V. Let 1Z° be the subset consisting of the elements a whose direction 
5 contains 0. 
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ei + e 2 / 













Figure 8. Decomposition for various generic points 

We now differentiate Req(<£>) with respect to (p. Using relations given 
in part (i) and (ii) of Lemma [9 .7\ we get 

^Req($)= ^($ / ,A,a,/3)-^^($/<0>,A/<0>,a/<0>,/3). 

We observe that the collection of a/ < > with a G 1Z° parametrizes 
all affine spaces admissible for $/ < <fi >. The collection TZ a g(^) may 
be larger than 7^ aff ($'), but, if a is not in 7^ a fj($'), then the contribution 
A(&, A, a, (3) is equal to 0. 

Hence, we obtain by induction that d^,(Leq(^>) — Req($)) = for 
any 0. 

Thus Leq($) — Req($) is constant. But by construction Leq($) — 
Req(<3?) is equal to zero on r, therefore the constant is zero. □ 

Example 9.8. We will give a decomposition formula for the system in 
Example 12.31 We recall the data: A = Zei © Ze2, $ = [ei, e± + e?\. 

We will compute i3($,A)(t>) using the decomposition formula for v 
in the tope defined by the inequalities 0<U2<l,l<Ui<2 and 
Vx — V2 > 1 (see figure[8]). 

We aim to demonstrate the dependence of the summands in the 
decomposition formula in Theorem 19.31 to the chosen generic point in 
a tope, though the value of £?(<&, A) (v) is clearly independent of this 
choice. We will thus decompose B(<&,A)(v) in two different ways, for 
two different choices of generic points lying in the same tope. 

For the first (resp. second) computation we choose a generic (3 = 
b\t\ + b 2 e2 (resp. 0') in the tope defined by < b\ < 1, < hi < 1 and 
b\ > b 2 , and further satisfying b\ + b 2 < 1 (resp. b\ + b 2 > 1). 
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Figure [HJ depicts two such choice of generic elements. We denote the 
projection of (3 and (3 1 to — e 2 © M(ei + e 2 ) by (3 and (3' respectively. 
Since these projections lie in different topes for the reduced system 
($ H s, A n s) their corresponding contribution to the sum in Theorem 
19.31 will be different. In fact, choosing as the generic point will en- 
force a nonzero contribution of the lattice point (1,0) in evaluating the 
distribution at a point v as depicted in Figure [HJ 

Computation with generic point f3: 

A)(v) = A($, A, V, P)(v) + A($, A, Re 2 + e h 0)(v) 

+A($, A, E( ei + e 2 ) - e 2 , f3)(v) + A($, A, {(1, 0)}, (3){v 

We now compute each summand using the formula in Lemma 17.31 
A($,A,Me 2 + e 1 ,f3)(v) = Ber(e 2 , t(/3 )) * T({e x , e x + e 2 }, (5\){v) 

= 5 {m * Res, =0 (((1/2 - dx 2 ) ■ ( -^gt^_ ) J 

= l(vi-l){-2v 2 +v 1 ) 



^($,A,R( ei + e 2 ) -e 2 ,p){v) = Ber(ei + e 2 , r(p )) * T{{e x , e 2 }, fii)(v) 

= 5 ( o,-i) * -Res z=0 (((1/2 - dx x ) ■ C-^-z^+z) )^ 
= -\(-Vi - v 2 )(vi -v 2 -l) 

A($, A, {(1, 0)}, 0)(v) = 5 (m * -T(e x , -e 2 , ei + e 2 ) = -fa - 1 - v 2 ) 
Using the computation in example 12.31 for A(<&, A, V,/3)(t>), we get 

jB($, A)(u) = -Uv 1 -2v 2 ){v x -l + v 2 ){2v l -l-v 2 ) + \{v 1 -l){-2v 2 + v 1 ) 

+ 2 (V1 + ^)fa - V 2 - 1) - fa - 1 - V 2 ) 

= -§fa - 1 - 2« 2 )(2« 1 - 3 - u 2 )fa - 2 + v 2 ) 

Computation with generic point (3': 

B{<$>,A){v) = A($,A,V,p r )(v)+A{$,A,Re 2 + e 1 ,P r )(v)+ 
A{^,A,R{e 1 + e 2 )+e 1 ,/3'){v). 

The first two summands in the decomposition above are already com- 
puted. The third summand equals: 

A(®, A, Rfa + e 2 ) + e u 0')(v) = --(2 -v x - v 2 )(v x - 1 - v 2 ). 
We then have 

B($,A)(v) = -Uv 1 -2v 2 ){v 1 -l + v 2 ){2v 1 -l-v 2 ) + l{v 1 -l){-2v 2 + v 1 ) 
-|(2 - vx - v 2 )(vx - 1 - v 2 ) 
= -§fa - 1 - 2^X2^ - 3 - u 2 )fa - 2 + v 2 ) 

as expected. 
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In the affine case, we define 

A($, A, a, r, u) = Ber(l> n s, r) * T($ \ s, u). 
We have a decomposition formula analogous to Theorem 19.31 
Theorem 9.9. Choose 6 V sufficiently generic. Then we have 

S(<£,A)=^(<£,A,a,/?). 

a 

The proof is precisely in the same line of arguments with that of 
Theorem 19.31 
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